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Abstract 

Glasses are known to exhibit quantitative universalities at low temperatures, the most striking of which is 
the ultrasonic attenuation coefficient In this work we develop a theory of coupled generic blocks with 

a certain randomness property to show that universality emerges essentially due to the interactions between 
elastic blocks, regardless of their microscopic nature. 
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1. Introduction 

Starting with the pioneering study of Zeller and 
Pohl experimental work over the last 40 years 
has shown conclusively that the thermal, acoustic, 
and dielectric properties of virtually all amorphous 
materials are not only qualitatively different than 
those of crystals, but also show a trulyamazing de- 
gree of quantitative universality 0, H, HI]. The the- 
oretical interpretation of the low temperature data 
on amorphous materials has for the last four decades 
been dominated by the "tunneling two state system" 
(TTLS) model. This model gives an attractive qual- 
itative explanation of the characteristic nonlinear ef- 
fects observed in ultrasonic and dielectric absorption 
(saturation, echoes, hole burning). In addition it 
predicts a frequency and temperature dependence of 
the ultrasonic absorption Q _1 (oj,T) (here defined as 
1~ 1 X/2tt 2 in terms of the phonon mean free path and 
wavelength) which appears to be in fairly good agree- 
ment with the experimental data; in the present con- 
text we note that Q^ 1 is predicted to be independent 
of u) and T, and the same within a factor of 2 in 
two regions of the parameter space, in which, at least 
in terms of the model, the physics is very different, 



namely the high frequency "resonance regime" (fc^T, 
ft )t <C io, where r is a characteristic relaxation time 
of the thermally excited TTLS's) and the low fre- 
quency "relaxation regime" (w <C h/r <C fcsT). This 
prediction appears to agree reasonably well with the 
data (see figures 2 and 3 in Q). Finally, at the cost 
of introducing a fairly large number of fitting param- 
eters, the model can reproduce most of the experi- 
mental data in the low-temperature regime reason- 
ably quantitatively. 

Nevertheless, there are a number of problems with 
the TTLS model. First, while in a few cases (such 
as KBr-KCN solutions 5]) it is possible to make a 
plausible identification of the "two level systems", 
in most amorphous materials their nature remains 
a matter of conjecture. Secondly, the model as such 
says nothing about the behavior at intermediate tem- 
peratures (1K-30K) which also shows a very strong 
degree of qualitative (though not quantitative) uni- 
versality. A third difficulty relates to the striking 
quantitative universality and small numerical value 
of the quantity (3 -1 (w); whether observed directly or 
inferred (in the "resonance regime") from the coeffi- 
cient of the log(T/Xo) term in the ultrasound velocity 
shift, this quantity has the value (3 ± 2) x 10~ 4 for 
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almost all non-metallic glasses measured to date 
While the TTLS model contains enough independent 
fitting parameters to "explain" this numerical result, 
the explanation requires a degree of statistical coinci- 
dence between these parameters which has no obvious 
basis in the model, and is prima facie nothing short 
of amazing. Finally, the model in its original form 
neglects the fact, which is emphasized below, that as 
a result of interaction with the strain (phonon) field, 
the TTLS's must acquire a mutual interaction |(|; 
while there exist theoretical approaches which take 
this feature into account and even use it 7] to at- 
tempt to account for the small universal value of 
it is not obvious that at this end of the process the 
TTLS structure is preserved, so that a question of 
self-consistency may arise. 

In Q, Q the conjecture was made that if we start 
from a very generic model in which at short length 
scales there is a nonzero contribution to the stress 
tensor from some non-phononic degrees of freedom 
whose only necessary feature is that their spectrum 
is not harmonic-oscillator-like, and take into account 
their phonon-mediated mutual interaction, we will 
recover at long length scales a picture which repro- 
duces most, if not all, features of the experimental 
data below 1 K. The goal of the present paper is 
quite modest: To attempt a somewhat more quanti- 
tative justification of this conjecture with respect to 
one specific feature, namely the (near) frequency in- 



dependence and small universal value of Q in the (1= 
regime ksT, H/t -C u> (i.e. what in the TTLS model 
is known as the high frequency resonance regime). 
We do not attempt to discuss here the behavior of 
Q^ 1 in other regimes (except for H/t <C uj < &bT), 
non-linear effects or (except briefly at the end of sec- 
tion 2.2) the intermediate-temperature behavior. 

We believe that it is one of the strengths of the 
present work that our result does not rely on ad- 
justable parameters, or the existence of other micro- 
scopic (unmeasurable) universal ratios to explain the 
observable one [HElHHEl (though cf.[7] ) . The 
only two inputs on which our outcome depends sen- 
sitively are the ratios ci/ct and Xl/Xt (cf- below for 
details of the notation) both of which are observed 
experimentally to vary little between different amor- 
phous systems (also cf. Appendix- A). Our third in- 



put ro, which is the size of a "microscopic amorphous 
block" (defined below) only enters into our equations 
logarithmically. 

The organization of the paper is as follows: In sec- 
tion 2 we define our model and introduce the central 
object of our study, namely the dimensionless stress- 
stress correlation, whose thermally-averaged imagi- 
nary part is the measured ultrasonic absorption Q~ x . 
In section 3 we carry out a real-space renormaliza- 
tion calculation of the average of Q'^w) over the 
frequency to and the starting state m (for details of 
the notation see below) and show that it vanishes 
logarithmically with the volume of the system and, 
for experimentally realistic volumes, has a surpris- 
ingly small value, ~ 0.015. In section 4 we argue on 
the basis of a calculation up to second order in the 
phonon induced interaction, that the functional form 
of Q _1 (a;) at T = should be (lnw)" 1 , and that 
when we combine this result with that of section 3, 
the numerical value of Q^ 1 for experimentally rele- 
vant frequencies should be universal up to logarith- 
mic accuracy and numerically close to the observed 
value 3 x 10 -4 . In section 5 we assess the extent 
to which our calculations are consistent with experi- 
ments in the (linear) resonance regime. In section 6 
we attempt to assess the significance of our results. 

Throughout this paper we set H = ks = 1- The 
notation a denotes a "typical" atomic length scale. 
The suffix a = l,t denotes the phonon polarization 
(Z=longitudinal, i=transverse) . 



2. Formulation of the Problem 

Consider a cube of an arbitrary isotropic amor- 
phous material, with side L which is assumed large 
compared to "microscopic" lengths a such as the typ- 
ical interatomic distance, but is otherwise arbitrary. 
We define for such a block the strain tensor in the 
standard way: If u(r) denotes the displacement rela- 
tive to some arbitrary reference frame of the matter 
at point r, then 



dx4 



dUj 

dxi 



(1) 



(Note that the anti-symmetric part of the tensor 
dui/dxj corresponds to a local rotation; since a spa- 
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tially uniform rotation costs no energy, any terms in 
the Hamiltonian associated with this part will be of 
order higher than zeroth in the spatial gradients, and 
hence for the purposes of the ensuing argument irrele- 
vant in the renormalization-group sense; we therefore 
neglect any such terms in the following). 

Let us expand the Hamiltonian of the block in a 
Taylor series in the strain e^ : 

H = H +J2e i jf ij +O(e 2 ) (2) 

ij 

where the stress tensor is defined by 

= dH/deij (3) 

Note that in general, in a representation in which 
Hq is diagonal, will have both diagonal and off- 
diagonal elements. 

As usual, we can define the static elasticity modu- 
lus X > a fourth order tensor, by 

xf M ^V-'(d{f l] )/de lj ) eq 

= V- l {d 2 H/de l3 de kl ) eq (4) 

where V = L 3 is the volume of the block and the 
suffix "eq" denotes that the derivative is taken in the 
thermal equilibrium state (both sides of eqn© are 
implicitly functions of temperature T). Since by defi- 
nition the properties of an isotropic amorphous mate- 
rial must be invariant under overall rotation, symme- 
try considerations constrain Xij-ki t° have the generic 
form 

XijM = (xi - 2xt)Si 3 Ski + Xt(8ikSjl + SuSjk) (5) 

where xi an d Xt ar e the standard longitudinal and 
shear elastic constants; in the approximation of an 
elastic continuum, these are related to the veloci- 
ties c; and Ct of the corresponding longitudinal and 
transverse sound waves (of wavelength A such that 
a«A«I) bj0 

Xi.t = pelt (6) 

1 Note that despite the notation \ has the characteristics 
of a "stiffness" (~ inverse susceptibility) rather than a "sus- 
ceptibility" 



where p is the mass density of the material. Such an 
approximation however throws away all the effects of 
interest to us, as we shall now see. 

2.1. The Stress- Stress Correlation Function 

Let us separate out from the Hamiltonian, the 
purely elastic contribution H e i, namely, 

Hei(eij) = const. + J ^d 3 r^xfj M e io^f) e u{r) 

ijkl 

i 

(where it is understood that the velocity is slowly 
varying over distances a, as above). Similarly we de- 
fine the "elastic" contribution to the stress tensor Ty 

by 

i&^ExgUi (8) 

ijkl 

(In above, (and ui) should strictly be treated as 
operators, but we prefer not to complicate the nota- 
tion unnecessarily). Then quite generally, we have 

H{e ij ) = H a (e ij ) + H'(e ij ) (9) 

where for the moment the "non-phonon" term 
H'feij) is completely general (in particular we do not 
assume it is necessarily "small" compared to H e [ ) . In 
analogy to @ and © we can define a "non-phonon" 
contribution to the stress tensor by 

H' =H^+Y / e ij f( j + 0(e 2 ) (10) 

ij 

% = dH'/dea (11) 

From now on we shall always take the quantities Hq 
and Tij to refer to the non-phonon contributions and 
accordingly omit the primes. Note carefully that the 
origin of in (|10p is not specified, and in particular 
it may have contributions from the phonon field. 

We can now define the quantity which will be the 
central object of our study in this paper, namely the 
(non-phonon) stress-stress correlation function (lin- 
ear response function) at different scales L. Consider 
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an externally imposed infinitesimal sinusoidal strain 
field, 

eij(r,t) = eijexp{i(q.r — cut)} + c.c (12) 

(eij real). This will give rise to a corresponding re- 
sponse of (Tij): 

(T y -)(r,t) = (T tf )exp{*(£?- wt)} + c.c (13) 

where (T^) is in general complex. Then we can de- 
fine the complex response function Xij--ki{q,u) in the 
standard wajH 

Xij>H ^ w ) = ^^{q,u) (14) 

We will usually omit the explicit g-dependence of 
X; it should be remembered that in general q and 
U)/ci t t are independent variables, cf.footnote ^ below. 

In practice we shall usually be interested in values 
of \q\ which are close to u)/ci t t and will therefore usu- 
ally omit the explicit q dependence of x (However cf. 
footnote El below) . 

Since for the purposes of the argument below, we 
shall be interested, at a geometrical scale L, in values 
of |g| which are of the order of L , it is not immedi- 
ately obvious that the symmetries of Xij-.M (q, w) allow 
its representation in the simple form analogous to ([5]); 
however, since it is clear that any complications as- 
sociated with this consideration are sensitive to our 
arbitrary choice of block shape, we will assume that 
a more rigorous (q-space) calculation will get rid of 
them, and thus assume that Xijkl(w) can indeed be 
represented in the form ([5]), thereby defining "longi- 
tudinal" and "transverse" response functions xz,t(w). 

All the above considerations are independent of the 
scale L of the block considered, provided only that 
this is large compared to atomic scales a. Let us now 
for a moment specialize to values of L of the order 
of the wavelength of the phonons studied directly in 
glasses. In view of the small values of Q" 1 and related 
quantities observed experimentally in this regime, it 



2 For brevity we omit some technical complications involved 
in the precise definition of \\ these are very standard, see e.g. 

m 



is very plausible to assume that the coupling between 
the phonon and non-phonon degrees of freedom (part 
of the second term in (ITU1) ') is a "weak" perturba- 
tion on the phonon dynamics as described by H e i. 
With this assumption it is straightforward to calcu- 
late the dimensionless ultrasonic attenuation Q~ 1 (w) 
of a phonon of frequency uj. Omitting the details of 
the derivation, we just quote the result 

Q-\uj) = A/-V27T 2 = Imxa(w)/(7rp<&) (15) 

where the 4 th rank tensor quantity Imx(w) (which we 
shall need in full generality below) is given explicitly, 
in the representation in which Hq is diagonal, by the 
formula 

lm Xlj M (w) = ^2 PmX&ll M ( 16 ) 

m 
n 

x S(E n - E m - u) (17) 

where m, n denote exact many-body eigenstates of 
Hq, with energies E m , E n , and p m is the probability 
of occurrence of initial state m (in thermal equilib- 
rium we of course havep m = Z" 1 exp{— f3E m }, where 
Z is the partition function). 

Since the above formulation is very generic, it is 
clear that the standard TTLS model must be a spe- 
cial case of it, specified by particular choices of the 
matrix elements of Hq and T. Without investigating 
these choices in detail (cf. below), let us note that 
according to this model the form of <5~ 1 (a;) is 

Q~V) = g,; / 1 tanh(w/2T). (18) 

This form seems to agree reasonably well with ex- 
periment, with Q^j ~ 3 x 10~ 4 . Note however 
that if (fT8|) is valid for all T, then in view of the 
Kramers-Kronig formula relating Q^j to xoq , the lat- 
ter quantity actually diverges logarithmically in the 
limit T — > 0. This might suggest (though it does not 
of course prove) that the true T = form of Q~ l is 
actually a weakly decreasing function of decreasing 
w; we will provide some evidence for this conjecture 
in section 4. 
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2.2. Phonon- Induced Stress- Stress Interaction and 
Real-Space Renormalization Program 
Consider a whole set of blocks, each of size L and 
described by the generic Hamiltonian (ITU1) ; for the 
moment we work at zero temperature. Since the 
quantity contains a contribution from the phonon 
field, it is clear that exchange of virtual phonons will 
give rise to an effective ("RKKY"-type) coupling be- 
tween the stress tensors of two different blocks 1 and 
2, and since the phonon system by itself is harmonic, 
this will have the generic form 

H^ = ^J drj dr"^2^Mr-r)T l3 (r)T kl (r'). 

ijkl 

(19) 

The function A^^r — r*) is calculated in the paper of 
Joffrin and LevelutQ; while the paper is formulated 
explicitly in the language of the TTLS model, it is 
clear that their form is generic. In the limit of large 
\f— f f \ (see below for the meaning of "large") it has 
the following form (l6| . where n is the unit vector 
along f — f 1 , 

(20) 

Ajjfej = —(Sji - 3rijni)Sik + (1 - Cf/cf)[~5 tJ S k i 
-SikSji - SuSjk + 3(Sijnkni + 5i k njni + SuUkUj 
+S jk nini + S^UiUk + 6 kl mnj) - I5nin k n k ni] 

(21) 

An approximation which should be adequate for our 
purposes is to replace r — f in ([T9| by R\ — R2 when 
R s denotes the center of block s, and that J v Tij{r)dr 
is the uniform stress tensor of the block. Then the 
total Hamiltonian of the N coupled blocks is, in an 
obvious notation, 

N 1 N 

(22) 

Eqn ([22|) then represents the Hamiltonian Ho of the 
"super block" (of side ~ N^-^L) composed by the iV 



blocks of side L; we can then define the stress ten- 
sor Tij for this super block and iterate the procedure 
until we reach the experimental scale. Note that be- 
cause of the factor proportional to L~ 3 in A«j.; and 
the fact that the correlation functions of the are 
defined intensively to have the same factor of L~ 3 , 
the procedure is scale invariant (in 3 space dimen- 
sions) and we might hope that it will iterate to a 
fixed point at large length scales. 

The program we would ideally like to implement, 
therefore, is to start with given forms of Hq and 
at some "short" length scale, introduce interaction 
between the corresponding "small" blocks according 
to (|22p . diagonalize the Hamiltonian Hn, obtain the 
corresponding stress tensor dHj^/de\j and iterate 
the procedure up to the experimental length scales. 
In particular, we would be interested in the extent 
to which the "output" forms of and at the 
experimental scale, are independent (or not!) of the 
"input" forms at the starting length scale. In prac- 
tice, a meaningful implementation of this problem 
appears to require massive computational resources, 
so in the rest of this paper we will concentrate on a 
few results which can be plausibly obtained by ana- 
lytical techniques. 

First, however, we need to discuss the question of 
an appropriate choice of L for the "input" blocks. 
It is well known that the temperature dependence 
of the thermal conductivity of amorphous materials 
changes at a temperature of the order of IK; in fact, 
most such materials show a pronounced "plateau" ex- 
tending very roughly between 1 and 10K. At 10K the 
dominant phonons have wavelengths of the order of 
50A, and in [§] it is argued that this is just the scale 
at which we get a crossover from "Ising" to "Heisen- 
berg" behavior (formally, at smaller scales the ap- 
proximations used in obtaining the simple R~ 3 form 
of A, cqn (|2"01) breaks down). Thus, we tentatively 
take the "starting" block size, which we will denote 
ro, to be ~ 50A(which is still comfortably greater 
than a). Since the results obtained in sections 3 and 
4 depend only logarithmically on ro , they will not be 
particularly sensitive to this choice. 
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3. Average value of Q 1 

To introduce this section let us consider the follow- 
ing heuristic argument (cf. Imagine a set of 
spin-like objects (in three spatial dimensions) which 
do not have direct mutual interactions but do inter- 
act with the phonon held via a term proportional to 
the strain, with some coefficient 7. As a result of 
this interaction, we will get an effective "spin-spin" 
interaction which is schematically of the form g/r 3 
times some angular factor which is of no great inter- 
est in this context, where g = ri^ 2 /(pc 2 ) with r\ a 
dimensionless number of order 1 (cf. below). 

Assuming that for the resulting effective Hamil- 
tonian it is possible to define some kind of single- 
excitation density of states P which is to a first ap- 
proximation independent of E, it then follows on 
dimensional grounds that P must be of the form 
const. Now in this model, the dimensionless ul- 
trasonic absorption Q" 1 is simply (7r/2)(7 2 / pc 2 )P, so 
it follows that Q _1 is "universal" (i.e. independent 
of 7, p and c). Moreover, since all phonon "colors" 
(modes) contribute to g, while only one color is ab- 
sorbed at a time, it follows that Q _1 ~ rf^ 1 ~ rT 1 
when n is the number of phonon "colors" . 

What we would like to do in this section is to try to 
(a) generalize this argument to a more generic model 
which does not necessarily assume "single-particle" 
excitations (with or without a constant density of 
states) (b) take into account quantitatively the ex- 
istence not only of different phonon "colors" but of 
different stress tensor components, and (c) argue that 
steps (a) and (b) alone lead to a surprisingly small 
value of Q^ 1 (though not small enough to explain by 
themselves the experimental data) 

To this end, it is convenient to define for a block of 
size L quantities Xa( u ) m terms of Xijkii 1 ^) ( a — 
l,t) in a way exactly analogous to that done in 
and further define an average Xa of X™( w ) over both 
frequency u and initial state m, and an associated 



3 It may be objected that the individual XiJfci( aJ ) do not 
necessarily have this symmetry; however, since all physically 
relevant quantities involve averages over m, we shall neglect 
this complication. 



dimensionless quantity Q a 1 as follows: 

fU 

Xa = N-'U, 1 £ / ^xi m V - E m ) (23) 
m Jo 

Oa 1 = {Kpcl)- l X a , (24) 

where N s is the number of levels below the cutoff Uq. 
The choice of Uq and the range of the sum over m 
will be specified below. 

We now consider two statistically identical blocks 
1 and 2, and write the complete Hamiltonian H of 
the coupled system in the form 

H = Hq + V (25) 

with Ho and V defined by the relevant special case 
of eqn (|2"2"]) . Since by hypothesis the matrix elements 
of Hq^ and 7^ (and thus of H and V) are uncor- 
rected, we may write 

TWT 2 - Tri? 2 = Trf 2 (26) 

provided, of course, that the three traces are over the 
same manifold of the complete many-body Hilbert 
space. The question of the optimal choice of this 
manifold is a delicate one. We know that for energy 
scales E n>m > Uq = hc/L and length scales \R\ < L, 
the assumptions that give rise to the simple form in 
cqn (|2"0j) no longer hold, and (m\V(R)\n) becomes a 
rapidly oscillating function of m,n, \R\. We handle 
these oscillatory effects by assuming that the interac- 
tions have no effect on E n > Uq, and the correspond- 
ing \n). Thus we take the relevant manifold to be 
simply spanned by those (many-body) eigenstates of 
Hq ' whose energies Eq are less than Uq. By restrict- 
ing ourselves within a low energy manifold, we are 
implicitly assuming that the interaction matrix ele- 
ments that mix in high-energy states into the man- 
ifold are negligible. To be quite explicit, we assume 
that the matrix element V mn is negligible whenever 
one of the eigenenergies m, n of H is less than Uq 
and the other greater; call this assumption (which we 
will carry over to subsequent iterations of the renor- 
malization process) assumption (A) . It seems unlikely 
that this approximation will influence the low- energy 
predictions that follow. In addition, we will see that 
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regardless of the length scale, Uq drops out of the cal- 
culation (cf. eqn([35|). and we need not worry about 
its precise value or L dependence. We hope to re- 
lax assumption (A), and possibly others in a future 
study. 

On general grounds, we may assume that the den- 
sity of eigenvalues E n of the interacting and nonin- 
teracting systems have the generic many-body form 

oo 

p(E) = a nE n (27) 

with dominating powers much larger than unity (in 
general, the dominating exponent of an extensive 
many-body system must be proportional to the num- 
ber of particles). While the actual many-body den- 
sity of states we deduce from the specific heat data, 

p(E) = p exp{const.(7V£) 1/2 } (28) 

is consistent with this requirement, as we will see, 
none of our results will be sensitive to the precise 
choice of a n . Because of the high power in eqn(|27|). 
and taking into account that the number of states 
in the superblock 1 + 2 is N 2 , we can evaluate 
TrH 2 w N 2 U^ (cf. Appendix-B). Moreover, if H 2 
is diagonalized within the manifold, it will have some 
maximum eigenvalue U , and by the same reasoning 
we expect TrH 2 w N 2 U 2 . Thus, 

TrH 2 - TrH 2 = N 2 {U 2 - t/ 2 ) (29) 

What of the RHS of ([26])? Here, there is a tricky 
point. We can evaluate the trace either in the eigen- 
basis of Hq or that of H. In the first case, it is clear 
that TrV 2 is proportional to N 2 Uq (Qq ) (where 
we now add the suffix to indicate that the quantity 
in question refers to the original size-L blocks); omit- 
ting the details of the algebra (cf. Appendix-B), we 
simply quote the result 

Tr^ 2 = CKN 2 U% (Qo 1 ) 2 (30) 

where we have defined a geometrical factor C = 
V 2 /[16it 2 (Ri — R2) 6 ], which for two side-by-side 
cubes is 9/167T 2 , and where the factor K is given by 

K= (8/3)[-3 + Ap+16q(q + p + qp- I)] (31) 



which for typical (experimental) values of q = 1 — 
c t/ c i ~ O-^ (cf- Appendix-1 for a theoretical justifi- 
cation) and p = Imx;/Imx t m 2.6, is 122. We note 
that the Qq 1 in ([50)) actually refers to the transverse 
coefficient Q^ 1 , the suffix of which we drop for nota- 
tional simplicity. 

It is also possible to evaluate TrV 2 in the eigenba- 
sis of H, in terms of the "renormalized" x' s °f the 
individual blocks 1 and 2 defined by taking, in the 
definition (TIT)) , the energies E m and states |m) to be 
the eigenvalues and eigenstates of H rather than Ho 
(recall that thanks to assumption (A), H can be di- 
agonalized within our submanifold) . It is clear that 
all the geometrical factors, etc., are identical in the 
two cases, so we get 

W 2 = CKN 2 U 2 {Qr 1 ) 2 (32) 

where (Qr 1 ) is the averaged absorption of one of the 
two blocks. We now face the difficulty that this is not 
necessarily the averaged absorption of the superblock 
1+2 (or even related to it by a simple numerical fac- 
tor), because the definition of the latter involves the 
squared matrix elements of the total stress tensor of 
the superblock, {f^'+T^') and thus contains terms 
like HT^V)H^ 2 V) (where | m), \ri) now denote 
eigenstates of H); while such terms were originally 
(in the absence of V) uncorrelated, it is not obvi- 
ous that they remain uncorrelated after V is taken 
into account. We shall, however, argue that on av- 
erage those terms are likely to be small compared to 
terms of the form \{m\T^\n)\ 2 , because V involves 

all tensor components of and while the cor- 
relation only involves the same component of 
and T( 2 );call this assumption (B). If this argument is 
accepted, we can identify the Q~ x in (|3"2")l with the 
physical inverse absorption of the superblock, which 
we denote simply Q^ 1 . 

Thus, putting together eqns.O,® and ([52]). we 
have 

N 2 (U 2 -U 2 ) = CKN 2 U 2 (Q^) 2 

= CKN 2 U 2 (Q^ 1 ) 2 (33) 

from which we can express Q _1 (the inverse absorp- 
tion of the superblock 1+2) in terms of that, Q^ 1 of 
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the original blocks 1 or 2: 



[1/iQo 1 ) 2 + CK] 



1/2 



(34) 



We can now repeat this procedure, by combining the 
superblock 1+2 (which we recall is the two cubes 1 
and 2 side by side) with a similar superblock 3+4, 
and finally combining the square structure so formed 
with a similar one to form a cube of side 2L. The only 
difference with the operation carried out explicitly 
above lies in the geometrical factor C, which is 4 
times larger than that in (I3U1) for step 2 and 16 times 
larger for step E@. As a result, the Q" 1 of the cube 
of size 2L is related to that of the original size-L one 
by (now denoted Ql for clarity) 



Q 



(2L) 



1 



189 
16tt 2 



K 



-1/2 



1 



1 -1/2 



(35) 



with K Q ~ 150. This completes the first stage of the 
iteration: We see that even in a single step of the it- 
eration the effects of the phonon-induced stress-stress 
interaction is to strongly renormalize the average ul- 
trasonic absorption downwards. 

We now consider the effect of iterating the step 
which led to (|35|). by combining eight cubes of size 
2L to make one of size 4L; for convenience we keep 
the original definition of the low-energy manifold 
(E m < wo), though other choices are also possible. 
Since the only point at which loq actually enters the 
result is implicitly in the definition of the "average" 
in Q^ 1 the scale invariant nature of the problem im- 
plies that all considerations are exactly the same as 
at the first stage, and we simply recover (|35p with the 
replacement of L by 2L and 2L by AL. Continuing 
the iteration up to a spatial scale R, we find 



Q-\R) 



Q 2 



(L) 



K \og 2 (R/L) 



-1/2 



(36) 



with the constant Kq weakly material dependent and 
approximately 150. It is remarkable that this formula 
predicts that at T = the ultrasonic absorption com- 
pletely vanishes in the thermodynamic limit. While 
we know of no particular reason why this behavior 
must be unphysical, in practice we would guess that 
for any finite ultrasound wavelength A, R would be 
replaced by a quantity of order >Q For experimen- 
tally realistic values R and the choice of L ~ 50A(cf. 
end of section 2) the value of Q _1 is approximately 
O.OljQ. While Q" 1 is material independent, its value 
is rather large compared to the experimental value 
in the MHz-GHz range. We think this is due to the 
contribution to the average, of the rapid increase in 
Q^ 1 at higher frequencies, as manifest in the thermal 
conductivity data around 10K. To obtain the experi- 
mentally observed absorption in the MHz-GHz range 
we need further considerations, one of which we now 
explore in section 4. 

4. The Frequency-Dependence of xi^*)'- Sec- 
ond Order Perturbation Calculation 

In the last section we obtained an expression (1361) 
for the value of (5 _1 (w) (or equivalently x(w)) aver- 
aged over the specified (low-energy) initial states m 
and over a frequency range U , but we could deduce 
nothing about the frequency-dependence of x- ln the 
present section we shall attempt to deduce some con- 
clusions about the frequency-dependence on the ba- 
sis of the specific ansatz (|38[) concerning the form 
of the "input" x at scale L. The reader should be 
warned that the argument of this section is heuristic 
and somewhat unorthodox. 

We consider the consequences (at T = 0) for x(uj) 
on taking into account the interaction term in V in 
(|22[) up to second order in perturbation theory for 
two neighboring blocks 1 and 2. In general there will 



4 To be sure, the approximation made in going from l |19ll to 
i ]22[) is somewhat less reliable than when combining two cubes 
side by side, and different ways of treating the sum (integral) 
can alter C; here, as at several other points, it would seem 
that a more reliable estimate of the numerical factors would 
require an explicit k-space calculation. 



5 To avoid confusion, we shall emphasize that in general 
(e.g. in considering the average over ui in J36}) the inverse 
length scale q = should be regarded as a different variable 
than ui/c; it is only when evaluating the predictions for the 
experimentally measured value of Q~ 1 that we should take 
q = u/c 

B Or slightly larger if we take R ~ A 



8 



be changes in (eqn ([T7|) ) due to (a) modifica- 

tion of the many-body density of states (DOS) and 
(b) modification of the matrix elements of the stress 
tensors Tjj (which, obviously, we need to calculate for 
the superblock 1+2, not for blocks 1 and 2 individu- 
ally). An important assumption which we shall make 
in this section is that, at least for the purposes of ob- 
taining qualitatively correct results, effect (b) may be 
neglected^ so that all we need to do is to calculate the 
shift due to V in the energies E n of the superblock. 
The first order term vanishes because of the lack of 
correlation in the matrix elements t£- and T^: the 
second order term can be expressed using eqns (|17|) 
and ([55)1 (when A is the correction to the energy of 
a superblock state E n = E ni + E n2 = lj) 

/oo poo 
/ q ~;»>')q^>") 

x6(U - \lu' + oj"\)(-oj' - cu")- 1 duj , duj" (37) 

where C and K are the quantities defined above in 
connection with eqnpip. To evaluate this formula we 
need an input form for Qq ■ For this we shall make 
the simplest possible choice which is consistent with 
our general assumptions, namely the "random" form. 

X { : n \u:)/^pcl) = const.6(E m +oj) 

= Qa 1 e{E m +u J ) (38) 

Note that eqn ([55|) does not specify the matrix ele- 
ments of Tij completely (it may or may not be one 
of the standard "random matrix" forms); we shall 
however postulate that T and H are uncorrelated. 
It should be carefully noted that the TTLS form of 
Xij^li is n °t a special case of ecin (j3"5|) ; this may be 
seen by noting that the form of Q~ 1 {oj,T) given by 
the latter is approximate!;^ 



7 One of us has in fact calculated llg|| the modification of 
x(u) to second order in V; the expression is an extremely messy 
sum of a large number of terms with different signs, and (un- 
surprisingly) its quantitative evaluation requires assumptions 
on the matrix elements of Tj which go beyond l |38H 

8 Note that because of the high power law in cqn l|27|l . the 
factor B(E m + u)) can essentially be neglected in calculating 
Q —1 (lj,T) (However it is possible that a more refined analysis 
will lead to a modification of this factor). 



Q- 1 {uj,T)^Q 1 (l-exp{-puj)) (39) 

which is different from eqn (|18[) , though not qual- 
itatively so. What we regard as most important, 
however, is that the form of Xi^lii 1 ^) i s very dif- 
ferent from that which would obtain were all the 
non-phononic degrees of freedom harmonic oscilla- 
tors; in the latter case we would simply get the re- 
sult Q^ 1 ^,! 1 ) = ind. of T, which is qualitatively 
different from the "saturating" forms (|39|) and (fT8|) . 
Crudely speaking, the ansatz (j38)) describes a model 
intermediate between harmonic-oscillator model and 
the TTLS one, but in some intuitive sense closer to 
the latter. It is our hope that use of a possibly more 
realistic form of the calculations of section 4 

will not lead to qualitatively different results. 

The expression (|37j) is finite for all lo (provided the 
principal part is correctly taken in the integral) . The 
ratio R(ui) of the DOS of the interacting superblock 
to the noninteracting one is given by 

R(u) = {1 + dA/duj)- 1 

= (l-CKQ 2 ln(U Q /u))- 1 (40) 

and in the absence of matrix element renormaliza- 
tion, this should be the expression for Q~ 1 (oj)/Q~ q 1 . 
However it is clear that we cannot take this result 
very seriously, since the expression (|4U|) is evidently 
negative for small (and even quite large, since the 
quantity CKQq 1 is 3> 1) values of lu. 

What has gone wrong? Let us consider the fol- 
lowing heuristic argument: Suppose that we could 
somehow neglect all higher-order effects of V and thus 
treat the result (1571) as exact. Then a negative value 
of R(lo) simply means that the order of levels in the 
region of ui has been inverted (i.e. each pair of levels 
crossed). This still give us a finite level density, but 
it is now given relative to its original value by the 
modulus of the expression (I4U1) ! Of course, when we 
introduce the higher-order effects of V, we will find 
that, barring pathology, we do not get level crossing, 
but rather the familiar level-repulsion effect; however, 
if we pretend for the moment that the overall effect of 
these higher-order terms are "small" relative to the 
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second-order ones, the resulting DOS should not be 
much affected by the lack of crossing. Consequently, 
we claim that at least over the large regime of uj for 
which CK 2 Qq 2 \u(U /uj) > 1, we have 



Q~ x {w) = [CKQa 2 ]n(Uo/u) - 1]" 



(41) 



which implies that for small uj 

Q~\lo) = const. (In Uq/uj)- 1 (42) 

It is clear that we can extend this result to take into 
account all pairs i,j in (|22p . in this case by a di- 
rect "single-shot" calculation. The result is to sim- 
ply replace the factor CK in (|4Tj) . as in (|36|) . by 
K \og 2 (R/L): 



Q-Hu) = 5 ^ (43) 

K Q- 2 log 2 (R/L)\n(U /u;)-l 

ps const. (Ln(i7o/w)) if uj < uj (44) 

While this form seems at least qualitatively consis- 
tent with the experimental data (see section 5), it is 
not even approximately universal (since the constant 
is inversely proportional to Qq 1 ), and even given a 
cutoff at ui ~ Uq does not satisfy (f36|) . Moreover, 
it is not at all clear that taking into account higher- 
order terms in V will not change the form of Q^ 1 
qualitatively. 

At this point, we adopt the following explicitly 
heuristic tactic: We know directly from ultrasound 
experiments that for uj <C Uq the frequency depen- 
dence of Q _1 (cj) is at least approximately consistent 
with that given in (I43[) ; we further know from the 
thermal conductivity data that for uj > Uq, Q -1 must 
have a much larger value, in fact, of order unity. We 
thus postulate for the overall frequency-dependence 
of Q _1 (w) the ansatjf] 



Q(w) _1 (w) 



1 



Qo + A\n(U /uj) 



(45) 



where Qq is not necessarily the same as Qq but is 
generally of the same order, and may, like Qq, be 



9 One may perhaps object that here (and in II46H ) Uo should 
be replaced by the renormalized quantity U. However, as we 
shall see, the conclusions are independent of Uq. 



system-dependent. We then find the value of A from 
the requirement that (|45|) be consistent with (f36|) . i.e. 
that 9 



Uo 



duj 



o Q + AlnUo/uj 



[(Qq 1 ) 2 +KQ\og 2 (R/L)]- 



KQ^log^R/L)]- 



(46) 



We see by a change of variables that the LHS of 
(|4"6")l is independent of the cutoff Uq and of the form 
A- 1 F(A- 1 Q ) where 



F(0 



dx 



o 1 - £ logx 



-e«Ei(-0 (47) 
ln(£ 2 e 2 ^)/2 (48) 



The approximation is valid for ^ < 1, and e 27 = 
3.17 ... is the Euler-Macheroni constant. Thus, pro- 
vided A~ 1 Qo < 1, we have 



Q = A^Iw^Qq/A) 



(49) 



note that this result is only very weakly system- 
dependent. If we now put Qq « 1 and require consis- 
tency of HI with ([H, we find that A w 350. Now, 
setting in (|45|) uj/(2tt) to be of the order 1MHz (a 
typical experimental value), we find 

Q^ilMHz) = 2.7 x 10~ 4 

which is precisely the "typical" experimental value. 

Since the ansatz (|45[) is only a conjecture, it is in- 
teresting to consider the somewhat more general form 

Q-\uj) = (Ql /S +A\u(Uq/uj))- s 

and inquire about the sensitivity of the value of 
Q _1 (lMHz) inferred from requiring consistency with 
(|36[) . For a general s this requires numerical calcula- 
tion; some representative results are shown in Table- 
1. We see that while the value of Q~ x is apprecia- 
bly sensitive to the (presumably system independent) 
parameter s, it is only very weakly sensitive to the 
system-dependent quantity Qq; thus, perhaps unsur- 
prisingly in view of the remarks in section 3, the "uni- 
versality" of the ultrasonic absorption is more general 
than the specific ansatz (|4"5")) . 
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s = 0.5 


s = 0.7 


s = 0.9 


Qo = 0.1 


0.0024 


0.0011 


0.0006 


Oo = l 


0.0024 


0.0010 


0.0004 


Oo = 10 


0.0024 


0.0010 


0.0003 



Table 1: The dependence of the value of Q(ui 
1MHz) to its functional form 



5. Discussion 

In section 3 and 4 we have attempted to draw some 
conclusions concerning the absorption of ultrasound 
at T = 0, in the MHz-GHz frequency range. How- 
ever, no existing ultrasound experiments has ever 
probed the regime uj 3> T. Thus, in this section 
we shall compare our predictions with data on re- 
lated quantities such as the temperature-dependence 
of the ultrasound absorption and velocity, and the 
thermal conductivity. For the purposes of this dis- 
cussion we will assume (cf section 3) that the many 
body states with ui < U are typical of the states as a 
whole, i.e. that the (renormalized) quantity xtjkii 1 - ) 
is "on average" independent of m, except for the fac- 
tor 0(E m +uj); this greatly simplifies the predictions 
for the temperature-dependencies of various quanti- 
ties. 

5.1. Temperature- Dependence of the Ultrasound Ab- 
sorption and Velocity 
In the standard TTLS model the temperature- 
dependence of the absorption, expressed in terms of 
Q _1 (w) as above, is given by 



Q-\u,T) = Q^ 1 tanh(w/2T) 



(50) 



(where Q^f is predicted to be independent of lo). 
In the present model the predicted dependence is 
straightforwardly calculated by taking into account 
the possibility of transitions with stimulated emis- 
sion as well as absorption of a phonon, and using the 
assumption above concerning x^ m ^(oj), we find 



i(w)(l-exp(-w/T)) (51) 



where now Q h ^ = const./ hx{U/ui). In Figure 1 we 
show a comparison of the two formulae (|50|) and ([51]) 



with experimental data [19( which probes the regime 
u) > T. Since most experiments, including that of 
ref. [19(, vary T rather than w, we note that when 
inferring the quantity Q~^j (i.e. the value of Q _1 (w) 
calculated in section 4) from the raw data, which is 
actually taken in the regime w < T we should divide 
the value of Qh / inferred within the TTLS model by 
a factor of 2. 

We next turn to the temperature dependence of the 
ultrasound velocity, Ac(T), expressed as a fraction of 
the velocity cq at some reference temperature ?cL3 
From the standard Kramers-Kronig relation for the 
stress-stress response function, we havj^l 



Ac 
Co 



p4 



= -A 



Q-\u,T), 



(52) 



Were we to insert in ^ the TTLS form Q~ 1 (lj,T), 
eqn (|50p . we would get the standard result for T oj 



Ac 

Co 



TTLS 



Q^HT/To) 



(53) 



Note that the TTLS requires a modification to its 
density of states n(E) — n a (l + aw 2 ) with free pa- 
rameter a to fit the data above ~ 0.5K (cf. [2(| and 
Figure 1). In the present model, were we to treat 
Q~hf as a constant, we would get, 



A(Q^/2)[- e "/ T Ei(- W /T) - e -/ T Ei( W /T)] 

(54) 



AC/C : 



which is precisely the same result as (|53[) in the 
u -C T limit. In the presence of [ln([/o/w)] _1 , the 
integration must be done numerically, and is shown 
in Figure 1 and compared with the TTLS result (with 
unmodified density of states) as well as experimental 
dated! [II. 



5.2. Thermal Conductivity 

It is well known [2l[ that the thermal conductivity 
in amorphous materials below IK is due entirely to 



10 Because of the very small value of Ac/c, the precise choice 
of To makes very little difference in the following. 

11 Some care is necessary with the sign; cf. the definition of 
X(^) in eqn.l |16l > 
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—x icr 4 




Figure 1: Normalized inverse mean free path 
(top) and velocity shift (bottom). The present 
theory (solid) is compared against TTLS without the 
uj 2 term (dashed) and experiment [ill [2(| (dots). 
While the functional forms predicted by both mod- 
els are qualitatively similar at low temperatures, the 
TTLS model must use an additional fitting function 
n(bj) = no(l + au! 2 ) for the density of states to resolve 
the discrepancy in fitting Ac/c data (bottom). 

phonons (any other degrees of freedom being effec- 
tively localized). If so, then it should be well approx- 
imated by the simple kinetic-theory formula 

K = ^CP h cl ph (55) 

where C% h is the phonon contribution to specific 
heat, c is a phonon velocity appropriately averaged 
over polarization {l,t) and l p h is a polarization and 
frequency-averaged phonon mean free path, which 
is proportional to Q. Since the frequency weight- 
ing heavily weights frequencies ~ 4fcsT, we may set 
Q(uj,T) approximately equal to its high frequency 
value and set l p h — const. Qhf(4T)/4T. Then for the 
TTLS model, l~ h L oc T, so the thermal conductivity K 
is predicted to be proportional to T 2 . In the present 



model, lph = Ixi(Uq/T)/T, so we predict 

K = const.T 2 ln([/ /T) (56) 

-a temperature dependence which is consistent with 
the usually quoted dependence of the experimentally 
measured thermal conductivity, namely K cx T 2_/3 , 
f3 ~ 0.05 — 0.2, and certainly fits it better than the 
TTLS prediction /3 = (see Fig.2). 



K(T)(Wcm"' K~') 




Figure 2: Temperature Dependence of Ther- 
mal Conductivity. The present theory (solid) 
K ~ T 2 In U/T is compared against TTLS prediction 
K ~ T 2 (dashed) and experiment [22| (dots) below 
the "plateau". 



6. Conclusion 

As already emphasized, our goal in this paper has 
been the relatively modest one of investigating how 
far it is possible to understand the linear ultrasonic 
properties of glasses in the "resonant" regime, and in 
particular the small and near-universal value of the 
zero-temperature dimensionless absorption Q _1 (oj), 
in terms of a model much more generic than the es- 
tablished TTLS model. 

We believe that we have been at least partially 
successful in doing so, in that we have shown that 
given assumptions (A) and (B) of section 3, a suit- 
able frequency average of Q{uj) is quasi- universal, in 
the sense that the only system-specific parameters on 
which it depends are the ratio q / c t and xi I Xt which 
fluctuate only by a factor of about 1.2 between dif- 
ferent materials (see Appendix, and refQ, Fig.l and 
3), and the high-energy cut-off Uq ~ hc/L, which en- 
ters only logarithmically, (b) that given the ansatz 
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(lnw) -1 for the frequency dependence of the 
absolute value of the absorption is close to the ex- 
perimental one. (c) that other properties related to 
(the temperature dependence of Q _1 and of 
the velocity shift, and the thermal conductivity) are 
consistent with our generic model. 

In this paper we have not attempted to discuss two 
other characteristic properties of glasses which are 
generally regarded as strong evidence for the TTLS 
model, namely the nonlinear behavior (saturation, 
echoes, hole-burning...) and the linear ultrasonic be- 
havior in the so-called "relaxation regime". We be- 
lieve that it is not impossible that a complete renor- 
malization calculation would produce, at the "out- 
put" stage, a form of x^ m \u) (suitably averaged over 
m) identical to that following from the TTLS model, 
in which case it would of course account for all the 
above phenomena equally well. 

Alternatively, if we assume that the output form of 
X^ m \(jj) is identical to the postulated input form (|35]l 
apart from the constant and a factor (lnw) -1 , it is in- 
tuitively obvious that at least some of the nonlinear 
behavior (e.g. acoustic saturation) will be qualita- 
tively similar to that in the TTLS model although 
the details will be different. The linear ultrasonic 
absorption in the low-frequency "relaxation" regime 
is however particularly puzzling, since the successful 
prediction of the experimentally observed frequency- 
and temperature-(in)dependence (and even the mag- 
nitude) by the TTLS model seems, at least at first 
sight, to be a consequence of parameter-dependencies 
very specific to that model. At the time of writing it 
is not clear to us whether it will be possible to repro- 
duce these dependencies in our more generic model. 
Of course, an alternative point of view would be that 
the TTLS spectrum is in fact always realized at the 
"input" stage (i.e. eqn ([3"8")l needs to be modified) and 
that its qualitative features are preserved through the 
renormalization process; the main role of this pro- 
cess, would then be to give rise to the universal small 
value of <5 _1 , and the final picture which would then 
emerge would be somewhat similar to that of Burin 
and Kagan Q. 

Finallywe note that the only property of amor- 
phous materials which we have invoked in this paper 
is the conjectured existence in them of an appreciable 



number of non-harmonic degrees of freedom. Thus, 
our arguments should apply equally well to "disor- 
dered crystals" (cf. [14]), with the slight modification 
that the angular dependence of h-ijki may be different 
in detail. However, we also note that nothing in our 
argument implies that all amorphous materials (or 
a fortiori disordered crystals) must have the "canon- 
ical" value of Q -1 ; if the input (scale-ro) value of 
Q _1 is appreciably smaller than this value, then so 
will be the output one! However, observation of an 
amorphous solid with the output Q~ x appreciably 
greater than the canonical value would tell against 
our hypothesis. 
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sity of Waterloo, and Harvard University for gener- 
ously providing accommodation during long visits. 

Appendix A. Transverse to Longitudinal 
Elastic Ratios 

The present model depends on the numerical values 
of two near-universal "inputs" , q /q and Im^t /Im^; . 
While these values are experimentally justified [3], 
we feel it may be instructive to add some theoretical 
justification at least to the trivial first ratio. 

Consider a many-body Hamiltonian that consists 
of an arbitrary inter-block interaction <fr(r) and ki- 
netic energy, which for small strains reduces to (J7J0 

Hel(eij) = const + L~ 3 / <f>(r)dr 3 (A.l) 



12 According to the virial theorem the expectation value of 
the potential energy is equal to that of the kinetic energy for 
a harmonic Hamiltonian for both longitudinal and transverse 
terms; thus considering the zero point kinetic energies does 
not change the ratio Xot/xot and will be omitted from our 
notation. 
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Further, suppose that the relative displacement Uij of 
two blocks i and j are proportional to the distance r*y 
between them. For example for purely longitudinal 
and transverse strains, we would have u x = e xy r y 
and u x = e xx r x respectively. Then, by definition, 



Xot,l 



L- 3 



(A.2) 



where i ^ j and i = j give \ot and xoi respectively. 
Letting d<j>/dr — due to stability, and taking the 
derivative d/deij using the distance proportionality 
assumption, we obtain \ot oc J r x r 2 dfl and xo; oc 
/ r^dCl, where f2 is the solid angle, and the constants 
of proportionality are identical. Since the speed of 
sound is related to the real part of the zero frequency 
response function according to ©, we get 



1 

71 



which is 6% larger than the experimental (average) 
value. 



Appendix B. Evaluation of the Traces 



if the terms for which n 3> 1 dominate the density 
of states (a premise required for the extensivity of 
many-body energy levels), we can use eqn (|B.l[) to 
write 

ETjn+l 
Cn sa U 2 N 2 (B.2) 
n + 3 s v ' 

n— 

The evaluation of TtHq is similar. Now let us turn 
to the trace of the square of 



V ab = 



3v ab pcf 



%*Hjkl J -ij 1 kl 



(B.3) 



where v a b — 47rr^ b /3 is the volume of a sphere defined 
by the inter-block separation r a b- 



M = TrV 2 = E A «« A 



i'j'k'V 



™i™2 ijkl 



(A.3) x (ni \f g \ mi )( ni \f? r \ mi ) (n 2 \f b kl \m 2 ) (n 2 |T fc % |m 2 ) 

(B.4) 



where \n s ) is the n th eigenstate of block s. The sum 
over states n s and m s is over the manifold described 
in the text, such that < E ms ,E nB < U. Thus, for 
all n s we may write the identity as 



In this appendix we outline the algebra involved 
in obtaining eqn (|2T))) and eqn (|3"2"j) (the derivation of 
eqn (|30)) is similar to the latter). 

Let us expand the non-interacting and interacting 
two-block density of states p(E) in powers of E, such 
that the the coefficient of the n th power is co« and 
c n respectively (cf. eqn([55])). For each set of coeffi- 
cients the normalization conditions require that the 
two body system has N 2 levels whether they interact 
or not 



v 



Jjno+l 

n "0 + 1 
n— 



V 

E< 

n=0 



u n+1 



= AT. 



(B.l) 



The trace in question is 



f 

TrH 2 = / LU 2 f{uj)dLU = V c„ 
J ° „ =0 



U 



n+3 



3' 



U-E„ 



S(E ria ~ E„ ls - uj')duj, 



insert it in the sum twice with s = 1 , 2 and evaluate 
the sums over ni,n 2 



M = 



^VaVb 



E 



V-E m , rU-E 



9(TTV ab pC 2 ) 2 J-E m J-E„ 

ijkl 
ijkl 

(B.5) 

= 9(lTVabpC 2 ) 2 U2N * E A iM A i'fkn'Xij-.i'fX b kl:k'l' 
i'j'k'V 

(B.6) 

Now let us decompose each factor of XijM into 
the only two independent components, as in eqnQ. 
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Since A^; only depends on q, and any Xij-.ki can be 
written in terms of Xt and p, it is obvious that the 
right hand side will be proportional to (Qt 1 ) ■ 

m = M!iwq-i) 2 (b.7) 

Vab 

where 

K = —2 ^2 ^ijkl^i'j'k'l'Xij.i'j'Xkl-.k'V (B.8) 

While the evaluation of the sum over all tensor 
components may seem complicated at first sight, the 
symmetries of Xiji'j'Xklk'l' simplify the problem con- 
siderably. Any term including a Xijkl with odd num- 
ber of equal indices (such as X1323 or X2223) van- 
ish. Furthermore, any term with indices {ijkli'j'k'l'} 
is equal to that with {i'j'klijk'V}, {ijk'l'ijkl} and 
{i'j'k'l'ijkl}. The sum is evaluated to yield K = 
(8/3) (-3 + 16(-1 + q)q + 4p(l + 2q) 2 ) « 122. 
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